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Abstract. We investigate here sums of triangular numbers f(x) := '^2,biT Xi where T n 

i 

is the n-th triangular number. We show that for a set of positive integers S there is 
a finite subset So such that / represents S if and only if / represents Sq. However, 
computationally determining Sq is ineffective for many choices of S. We give an explicit 
and efficient algorithm to determine the set So under certain Generalized Ricmann 
Hypotheses, and implement the algorithm to determine So when S is the set of all odd 
integers. 



1. Introduction 

In 1638 Fermat wrote that every number is a sum of at most three triangular numbers, 
four square numbers, and in general n polygonal numbers of order n. Here the triangular 
numbers are T x := £i^hll ) where we include x = for simplicity. The claim for four 
squares was shown by Lagrange in 1772, while Gauss famously wrote "Eureka, A + A + 
A = n" in his mathematical diary on July 10, 1796. 

Theorem (Gauss, 1796). Every positive integer is the sum of three triangular numbers. 

The first proof of the full assertion of Fermat was given by Cauchy in 1813. 

In 1917, Ramanujan extended the question about four squares to consider which 
choices of b = (bi, b 2 , b 3 , 6 4 ) satisfy b\x\ + b^x\ + b^x\ + b^x\ representing every positive 
integer. We shall refer to such forms as universal diagonal forms. He gives a list of 55 
possible choices of b which he then claims are the complete list of universal quarternary 
diagonal forms (54 forms actually turned out to be universal). 

In 1862 Liouville similarly proved the following generalization of Gauss's theorem. 

Theorem. Let a, b, c be positive integers with a < b < c. Then every n e N can be 
written as aT x + bT y + cT z if and only if (a, b, c) is one of the following: 

(1,1,1), (1,1,2), (1,1,4), (1,1,5), (1,2,2), (1,2,3), (1,2,4). 

In fact, the following simple condition determines whether a fixed set of positive 
integers (bi, . . . , bk) give rise to a sum of triangular numbers Y2i=i bkT Xk which represents 



every integer as shown by the author in (114l ). 
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Theorem 1.1. Fix the sequence b\ < ■ ■ ■ < bk G N. Then 

1. The sum of triangular numbers 

k 

f(x):= fb(x) := Y^ b i T ^ 
i=i 

represents every positive integer if and only if fb represents the integers 1, 2, 4, 
5, and 8. 

k 

2. The corresponding diagonal quadratic form Q{x) = ^2,bixj with Xi all odd rep- 

i=l 

k 

resents every integer of the form 8n + J2 b i with n > if and only if it represents 

i=i 

8 + 16 + $>, 32 + 40 + and 64 + 

i=l i=l i=l i=l i=l 

Recently, Conway and Schneeberger proved a very nice similar condition for positive 
definite quadratic forms whose corresponding matrix has integer entries, but without 
publishing their results. 

Theorem (Conway-Schneeberger). A positive definite quadratic form Q(x) = x t Ax 
where A is a positive symmetric matrix with integer coefficients represents every positive 
integer if and only if it represents the integers 1, 2, 3, 5, 6, 7, 10, 14, and 15. 

Bhargava gave an elegant simpler proof of the Conway-Schneeberger 15 theorem in 
(0), in addition to showing more generally that for any set S C N it is always sufficient 
to check whether Q represents a finite subset So, and showed the set So for the two sets 
S = {2n + 1 : n G Z + } and S = {p prime}. 

In this paper, we will consider a similar generalization of Theorem 11.11 

Theorem 1.2. Let a set S CN be given. Then there is a finite subset Sq of S such that 
f(x) represents S if and only if f represents Sq. 

A simple computer calculation leads us to conjecture a set So when S is the set of all 
odd integers, for example. 

Conjecture 1.3. A sum of triangular numbers f represents all odd integers if and only 
if it represents the integers 

1, 5, 7, 9, 11, 13, 17, 19, 25, 29, 35, 49, 89. 

Unlike in Bhargava's theorem however, current techniques are insufficient for compu- 
tationally determining a suitable So for most choices of S, due to ineffective bounds for 
the class numbers of imaginary quadratic fields. We shall briefly explain this complica- 
tion. Let f(x) = biT xl + b 2 T X2 + 63^3 be given such that / represents all of the integers 
in 5", but the corresponding (diagonal) quadratic form is not (spinor) genus 1. Then the 
corresponding weight 3/2 modular form with all odd can be written as an Eisenstein 
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series plus a cusp form. Siegel has shown that the Fourier coefficients (with bounded 
divisibility at the anisotropic primes) of the Eisenstein series grow like the class number 



(23). Siegel has also shown that the class number grows faster than n^~ e (22), but the 
bound was ineffective because Siegel showed the result by first assuming a certain Rie- 
mann hypothesis and showing the result, and then assuming the Riemann hypothesis 
was false, and getting a different constant, depending on the location of possible zeros. 



The best known effective bound is given by Oesterle (I16l ). but is only O(logn). After 
decomposing the cusp form into g\ + Q2 where g\ is in the space of lifts of one dimen- 
sional theta series and the Shimura lift of #2 is cuspidal, we note that the coefficients 
of <?2 grow slower than n5"i +£ by the work of Duke (0), and g\ is supported at finitely 
many square classes with the same growth as the Eisenstein series. Since the coefficients 
with bounded divisibility at the anisotropic primes of the Eisenstein series grow faster 
than the coefficients of every sufficiently large n with bounded divisibility by the 
anisotropic primes and outside of the support of the coefficients of gi must be repre- 
sented. However, we do not know the implied constant from Siegel's ineffective bound, 
so we cannot effectively determine when n is sufficiently large. 

Assuming the Generalized Riemann Hypothesis for Dirichlet L-functions, and using 
Duke's (effective) bound of rfi +e (0), we would have an algorithm to determine whether 
/ represents s or not. However, although Duke's result is effective, the author is un- 
aware of any paper where the constant is explicitly computed. Even assuming that the 
implied constant was 1, the bound obtained is entirely infeasible with current computer 



technology. Using an idea of Ono and Soundararajan (I18l ). and a generalization by the 



author (jl2l ). we will be able to determine an algorithm to determine the set So under the 
additional Generalized Riemann Hypothesis for L-functions of weight 2 newforms. For 
notational ease, we will refer to an integer which is locally represented at every prime 
as an eligible integer. 



Theorem (Ono-Soundararajan (Il8l)). Assume GRH for Dirichlet L-functions and GRH 
for L-functions of weight 2 newforms. Then the eligible integers not represented by 
Ramanujan's ternary quadratic form x 2 + y 2 + 10z 2 are precisely 

3, 7, 21, 31, 33, 43, 67, 79, 87, 133, 217, 219, 223, 253, 307, 391, 679, 2719. 

Using the ideas in ffl8l ). (fl2l ). and ffl3l ). we obtain the following result. 

Theorem 1.4. Assume GRH for Dirichlet L-functions and GRH for L-functions of 
weight 2 newforms. Let a set of positive integers S of nonzero density be given. Then 
there exists an explicit algorithm to determine the (unique) smallest (finite) set Sq such 
that a sum of triangular numbers f represents S if and only if f represents Sq. 

Moreover, this algorithm is computationally feasible for sets S which contain many 
small integers. 

We also implement the algorithm to prove Conjecture 11.31 under these GRH assump- 
tions. 
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Theorem 1.5. Assume GRH for Dirichlet L-functions and GRH for L-functions of 
weight 2 newforms. Then a sum of triangular numbers f represents all odd integers if 
and only if it represents the integers 

1, 5, 7, 9, 11, 13, 17, 19, 25, 29, 35, 49, 89. 

Moreover, this is the smallest such set with this property. 

Remark 1.6. We only need to assume GRH specifically for weight 2 newforms that occur 
in the decomposition of any f for which, under the GRH assumption we can show will 
generate all integers in S . However, to make our algorithm more efficient, we will make 
the GRH assumption in general so that we do not need to check very many coefficients 
of a large number of weight 2 modular forms (For more information on modular forms, 
please see ( fl"?b ). Thus, we note here that for determining the set Sq for S all odd integers, 
we could have avoided the GRH assumptions except in 3 cases at the added disadvantage 
of more computation. 

We conclude with a curious example. Consider the set 

S ■= {n G N : 2x 2 + 3y 2 + Az 2 = 8n + 9 has a solution}. 

Current techniques appear insufficient to determine the set S. However, we are able to 
prove the following. 

Theorem 1.7. The sum of triangular numbers f represents S if and only if f represents 

2,3,4,5,10,16,17,19,89. 

Notice that the above theorem is unconditional. As we will see in Section [31 GRH 
predicts that 

S = S := N\ {1,8,31}. 
However, without GRH we know of no algorithm to determine the set Sq. 

2. Existence of S 

We will begin by showing that fixing a subset S of the positive integers, that indeed 
checking a finite subset So will suffice. 

Proof of Theorem ll.ffl We will follow the basic argument of escalator lattices in Bhar- 
gava 's argument (J2|) for quadratic forms. Without loss of generality, we will denote 
the triangular form f = fb simply by the sequence of coefficients, b = [fox, . . . , bk], with 
b\ < ■ ■ ■ <bk- Note that we must represent the smallest integer G S, so it follows 
easily that foi < S0. For each choice of b±, we find the smallest integer sr^i G S which is 
not represented by [foi], and conclude that bi < sr^i. We recursively continue to build 
a tree of possible choices of bk depending on the previous choices of foi, ... , b^-i- Note 
that we never need the to choose the same integer fo« more than 3 times in one branch, 
since this then precisely represents every integer congruent to zero modulo fo, by Gauss's 
theorem. Whenever fo represents all odd integers, we will say that fo is a leaf of the tree, 
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since any arbitrary choice b' containing b as a subsequence will automatically represent 
every integer in S. For b G T not a leaf, we will denote by s b the smallest s 6 S not 
represented by b, and we will call s b the truant of 6. Thus, taking for our tree T 

S := {s b :be T} 

it follows easily that a triangular form / represents S if and only if it represents So by 
our construction of T. Moreover, such a choice of So is smallest possible and unique, 
since for the form b, s b is the smallest s G S not represented by b, noting that 

[bi, . . . , b k , s b + 1, s b + 1, s b + 1, s b + 2, s b + 2, s b + 2, . . . , (s 6 + l)(s& + 2) - 1] 

represents exactly every integer s G S other than s&, using Gauss's theorem that every 
integer is the sum of three triangular numbers. 

It therefore remains to show that So is finite. Since at each step there are only finitely 
many choices for bk, the breadth at each node of the tree is finite, so it suffices to show 
that the supremum of the depth is finite. To do so, we will consider all nodes at depth 
4. Since the breadth is finite, there are only finitely many such nodes, and only finitely 
many leaves of depth less than or equal to 4. Therefore, it suffices to fix one such node 
and show that the depth of the resulting subtree is finite. 

Let / be a triangular form of dimension at least 4, and consider Qodd to be the 
corresponding quadratic form with x\ odd. Then the n-th coefficient of / is the equal 

k 

to the 8n + ^fr^-th coefficient of Q Q dd- Therefore, the coefficients a/(n) are precisely 

i=l 

the coefficients of the modular form corresponding to Q dd- Hence, we may decompose 
these coefficients into coefficients of an Eisenstein series plus coefficients of a cusp form. 
Since the dimension of / is at least 4, the coefficients of the Eisenstein series grow faster 
than the coefficients of the cusp form as long as the Eisenstein series is non-zero. We 
note that there are finitely many congruence classes where the coefficients of both are 
zero, namely those integers not locally represented by the quadratic form. Hence, as 
long as the Eisenstein series is non-zero, there are only finitely many congruence classes 
and finitely many "sporadic" integers not represented by /, since the coefficients of the 
Eisenstein series are always positive. It is simple to show that the Eisenstein series must 



be non-zero, however, since by Siegel's local density formula (123l ). this Eisenstein series 
is given by the difference of the local densities with arbitrary Xi and the local densities 
with Xi even for some i. However, a quick check show that when p ^ 2 the local densities 

k 

are the same, and, since the integer 8n + Yl^k is locally represented modulo 8 with 

i=l 

all odd (namely %i — 1), the local density for arbitrary must be greater than when Xi 
is even except at finitely many congruence classes. 

Let one of the nodes of our tree, / G T be given of exactly dimension 4. Then 
there are only finitely many congruence classes and finitely many "sporadic" integers 
not represented by /, and hence at each step of our escalation, our next choice of s b G S 
must be one such integer. We may only escalate finitely many times for the "sporadic" 
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integers, and each time we escalate to include an element of one of the congruence 
classes, the congruence class is replaced with finitely many new "sporadic" integers at 
the next step. Therefore, since there are only finitely many congruence classes, we can 
only add finitely many new "sporadic" integers overall, and hence the subtree of / is of 
finite depth. □ 

3. Determining S 

We will describe the algorithm to determine the set Sq. For complete details of how 
to compute the bounds obtained given GRH, see 0) and (0). 



Proof of Theorem ~Lj From the proof in Section [21 it is clear that Sq is uniquely deter- 
mined by the tree T, so this algorithm is equivalent to determining the tree T, since it is 
a simple check to determine the smallest seS 1 which is not represented by a fixed form 
/. Constructing the tree as in the proof is also quite simple, so that the only remaining 
obstacle is determining whether a node of the tree is a leaf. Our task is thus equivalent 
to determining (effectively) which integers are represented by Q dd- If the dimension 
(and hence the depth in the tree) of / is at least 5, then, using the trivial bounds for 
the coefficients of the Eisenstein series and the cusp forms, we may effectively determine 
a bound beyond which every integer locally represented is globally represented as in 



Tartakowsky's work (1241 ). and local conditions are a simple check at the primes dividing 
the discriminant. For depth 4, the wonderful optimal bound for cusp forms of Deligne 
(jfl) and calculation of the anisotropic primes allows us to again effectively determine the 
set of integers represented by Q a dd (cf. Hanke (fiol)). For forms of depth 2 or less, we 
note that Qodd only represents a set of density zero, so that the form cannot be a leaf. 

It remains to determine whether a form of dimension 3 is a leaf or not. However, 
additional complications arise for ternary quadratic forms. First note that the inclu- 
sion/exclusion of theta series 

^Qodd : = 0Q(x,y,z)-8Q(2x,y,z) -&Q{x,2y,z) -&Q{x,y,2z) +^Q(2x,2y,z) +&Q{2x,y,2z) +9q(x,2 v ,2z) ~ 6 } Q(2x ,2y ,2z) 

has n-th coefficient non-zero if and only if Q dd represents n, where Q is the quadratic 
form without the odd restriction. We note that 9 := 0Q odd decomposes as follows: 

9 = (9 — 9 spin ) + (9 spin — 9 Gcn ) + (6> Gen ), 

where 9 sp i n is the inclusion/exclusion of the (Siegel) weighted average over the spinor 
genus and 9 te xtGen is the inclusion/exclusion of the weighted average over the genus. It 



is well known (cf. (J8|)), that 9q — 6> spin (Q) is the orthogonal complement of U, 9, 



spi 



(Q) 



^Gen(Q) £ U, and # Gcn (Q) is the Eisenstein series given by the local densities in f)23l). 

Firstly, as noted in the introduction, the coefficients of the Eisenstein series grow 
like the class number, so that we have made the assumption of GRH for Dirichlet L- 
functions. Additionally, we have assumed GRH for the L-series of weight 2 newforms. 
Additionally, there are what are referred to as anisotropic primes. That is, primes p for 
which the n-th coefficient of the Eisenstein series does not grow when n grows with high 
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divisibility by p. Luckily, a local condition allows one to check only finitely many primes 
(those dividing 2D, where D is the Discriminant) to determine these anisotropic primes. 
Additionally, the genus of the quadratic form may be broken into what are referred to 
as Spinor Genera. In each spinor genus, there are finitely many integers t, called spinor 
exceptions for which a subset of the square class tT? is not represented by the form. This 
comes from the fact that the form splits naturally into three parts, namely an Eisenstein 
series, a cusp form in the space spanned by lifts of one dimensional theta series (we will 
denote this space by U), and a cusp form in the orthogonal complement of U. The forms 
in the space spanned by the lifts of one dimensional theta series have coefficients which 
grow like in square classes tZ 2 , and hence with the same growth as the coefficients of 
the Eisenstein series. However, the coefficients of each lift is non-zero exactly at t times 
a square, and t is restricted by the level of the corresponding modular form. Schulze- 
Pillot has given an explicit algorithm to determine the full set of spinor exceptions for 
a given quadratic form (fl9l ). so this problem will be resolved by using Schulze-Pillot's 
algorithm. If the t-th coefficient of 6* sp i n is non-zero, then the growth of the coefficients 
of # S pin grows like the class number in this square class. Otherwise, by investigating the 
spinor norm mapping, Schulze-Pillot determines explicitly the integers m such that the 
tm 2 -th coefficient of 6 spm is zero and those for which the coefficient is equal to a positive 
constant times the Eisenstein series (cf. Schulze-Pillot fj20l)). for which we can reduce 
the problem to the case where t is not a spinor exception for the genus. 

Since we can determine the finitely many spinor exceptions, we may then ignore the 
part of the decomposition resulting from a cusp form in U, and it remains to give 
(effectively) a bound Ng such that the n-th coefficient of the Eisenstein series are larger 
than the n-th coefficient of the cusp from in U 1 - whenever n > Ng. 

To do so, we must have effective bounds for the coefficients of both the Eisenstein 
series and the cusp form. The Eisenstein series is 



n 

d 2 \n 



CgJ^h 



where Cg may be determined by the local densities in (j23l ). h(m) is the class number 
of the imaginary quadratic field Q(V— m) and a(4D,n) is a constant depending only on 



gcd(4D, n)where D is the discriminant (jlll ). 



For square free integers, we will use Dirichlet's class number formula (cf. (J5|) to rewrite 
the class number with the special value of the Dirichlet character X-n, 

h{ _ N) = V^L(X-N,l) 



For square free integers, a celebrated result of Waldspurger (125l ) allows us to rewrite the 



square of the absolute value of the coefficient of a Hecke eigenform as the special value of 



the (integral weight 2) Shimura lift (12 ll ) of the Hecke eigenform. Hence, we decompose 
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the cusp form into Hecke eigenforms and then use the Schwarz's Inequality, giving 

\ a O apin -0 G en( n )\ < 



where q is some explicitly computable constant given by a fixed n$ = n (mod 4.D 2 ), 
Gi is the Shimura lift, x 1S the Nebentypus (cf. (1171)) of the weight 3/2 cusp form, and 
L(Gi, is the special value of Gi twisted by ip. 

Hence, rearranging, if n is square-free such that the coefficient ao(n) is zero (and n is 
not a spinor exception), then 



i >A L(Gi,XX-n, 1 

cen* < > d— - 



^(X-nXap^j , l) 2 

where c# is an explicitly computable constant. Thus, it only remains to bound ^^^jr- 
n s effectively and with the implied constant given explicitly for some 5 < h. However, 



under the given GRH assumptions, an explicit bound is given in (1151 ) . and the details of 



the ensuing calculation for the implied constant are given in (ll3l ). 

For n not square free, the Hecke operators may be used to show an explicit bound 
for the squares part beyond which integers must be represented, away from the spinor 



exceptions. For more details, please see (1 121 ) or (1151 ) . Therefore, we can conclude that 
the set of square free integers not represented by this form may be exactly determined 
by checking up to the bound obtained, and hence we may (effectively and efficiently) 
determine whether this form is a leaf. For forms with small discriminant (say, less than 
300), the bound obtained is often well below 10 12 , and hence is well within the computing 
power of current technology. 

□ 

Remark 3.1. In practice, whenever a leaf exists at depth 3, we will determine in general 
the set of integers not represented by Q dd for the nodes at depth 3 (not just the leaves) 
and then note which of these integers remain at each step of the escalation, instead of 
using the arguments of Tartakowsky (24) at depth at least 5 or the bounds of Deligne (0) 
at depth 4- 



We now implement the above algorithm to show the set So given in Conjecture 11.31 
when S is the set of all odd integers, is the correct smallest such set under GRH. 

Proof of Theorem \1.5\ We will proceed by considering each node at depth 3, and deter- 
mining the corresponding subtree under these GRH assumptions. For notational ease, 
we will refer to the form 

r 

/ := 

i=i 

by [bi, . . . , b r ], and the corresponding quadratic form by (b\, . . . ,bk). 
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The forms [1,1,1], [1,1,2], [1,1,4], [1,1,5], [1,2,2], [1,2,3], and [1,2,4] represent every 
natural number by Liouville's theorem. Since [1, 1,3] is a genus one form, the integers 
not represented by [1,1,3] are precisely the integers n such that 



Therefore, we are missing the integer 17, and we must escalate to [1,1,3, &] for some 
k < 17. If n is not represented by [1, 1, 3, k], then it is not represented by [1,1,3], either, 
so 8to + 5 = 3 2r+1 (3£ + 2), so that n = 5 (mod 9) or n = 8 (mod 9), depending on 
whether r > or r = 0, respectively. But then, taking £4 = 1, it follows for n > k that 
n — k = 5 (mod 9) or n — k = 8 (mod 9), and hence it follows that 3 | k. If r > then 
it follows that k = 0, 6 (mod 9), and if r = it follows that k = 0, 3 (mod 9). 

For the case [1,1,3,3] we then note that the form Q = (1,1,3,3) is genus 1. Our 
inclusion/exclusion of theta series gives that n is represented if and only if the 2n + 2-th 
coefficient of 9q is positive, and it follows that every to is represented because the local 
conditions are always satisfied. 

For the cases [1, 1, 3, k] with = 6 or = 15 we have r > 0, so that we only need 
to consider to = 5 (mod 9), or in other words 8n + 5 = 3 2r+1 (3£ + 2) with r > 0. 
We check the cases n < 3k by hand and for n > 3k, the choice X4 = 2 shows that 
8 (to — 3k) + 5 = 3 2r +1 (3£' + 2), with r' > by congruence conditions modulo 9. Taking 
the difference and denoting R = min{r, r'}, we have 



But Vs(24k) = 2 and 3 2R+1 divides the right hand side, giving a contradiction since 



In the case [1, 1, 3, 12] we note that we have the truant 89 and r = from above. In 
this case we take x 4 = 1 to obtain to' = 8to + 5 = 3 2r ' +1 {3£' + 2) + 96 with r' > 0. We 
then have either r' — 1 and nf = 69 (mod 81) or r' > 1 and to' = 15 (mod 81). Assume 
r' > 1 and set X4 = 4 to obtain to' = 3 2r +1 (3£" + 2) + 960. Taking the difference, we get 

54 = 864 = 3 2r ' +1 (3£' + 2) - 3 2r " +l (3£" + 2) = -3 2r " +1 (3t" + 2) (mod 81). 

It follows immediately that r" = because otherwise the right hand side would be zero. 
But now we have to' = 3 3 (3£" + 2) + 960 = 42 (mod 81), which contradicts the fact that 
n! = 15 (mod 81). Hence we have only the case to' = 69 (mod 81) remaining. We now 
escalate to [1, 1, 3, 12, k] for 9 < k < 89 (although we have restricted our coefficients to 
be monotone, we include the case k = 9 for usage below). Since [1, 1,3, 12] represents 
every integer not congruent to 69 modulo 81, we are done when k ^ 81, and for k — 81 
the truant 89 remains. Furthermore, if we escalate further with k = 81 we will never 
obtain 89 since [k, k, k] represents precisely kN by Gauss's theorem. Furthermore, if we 
ever choose k = 81 and then escalate to another integer (such as [1, 1, 3, 12, 81, k}) we do 
not obtain a new truant because [1,1,3, 12, k] has no truants. Whenever this situation 



8n + 5 = 3 2r+1 (3£ + 2). 



24k = 8to + 5 - (8(to - 3*;) + 5) = 3 2r+1 (3£ + 2) - 3 2r ' +l {3£' + 2) 




R>0. 
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occurs henceforth we shall say that we are "stuck" at k = 81. Following the above, for 
[1, 1, 3, 9] we are stuck at k — 9. This concludes the subtree of [1, 1, 3], and since [1, 1] 
does not represent 5, we have included the entire subtree of [1, 1]. 

Our arguments for [1, 2, 5], [1, 2, 10], [1, 3, 4], and [1, 4, 6] will all be identical and similar 
to the cases above, so we will combine them together. We will demonstrate the argument 
for [1,2,5] and leave the other cases to the reader. First we note that the number of 
representations of n by [1, 2, 5] is the same as the number of representations of 8n + 8 
by (1, 2, 5) minus the number of representations of 8n + 8 by (4, 8, 20), since if any are 
even then all must be even, taking everything modulo 8. For simplicity, we will denote 
t[b](n) to be the number of times that n is represented by the sum of triangular numbers 
corresponding to b, tvm {n) for the number of times the quadratic form represents n, and 
r °b)( n ) f° r t ne number of times the quadratic form with all X{ odd represents n. So the 
above is simply 

£[1,2,5] (n) = r° (h2 5) (8n + 8) = r (li2 , 5 ) (8n + 8) - r (4,8,20) (8n + 8) . 

Now, r (4i8i 2o)(8n + 8) = r(i, 2 , 5 )(2n + 2). However, (1,2,5) is genus 1, so r^ 2 ,5)(rn) 



is 



given precisely by the local density (23). Checking the local densities, we see that 
r (1,2,5) (8n + 8) = 2r(i j2j5 )(2n + 2), since the local densities are clearly equal when p 7^ 2 
(taking the isomorphisms Xi 1— > 2~ 1 Xj), and for p = 2a simple computation shows exactly 
twice as many solutions modulo the same 2 power. Therefore, 

*[i,2,5](n) = r (1,2,5) (2n + 2). 

Again, noting that (1, 2, 5) is genus 1, we know that n is represented globally if and only 
if it is represented locally. The integers not represented locally by (1,2,5) are integers 
of the form 5 2r+1 (5n + m) where m is a non-square modulo 5. Therefore, it follows that 
if n is not represented by [1, 2, 5], then 5 | {n + 1). So [1, 2, 5, k] must be a leaf if 5 \ k. 
Since the truant for [1,2,5] is 19, we need only check k — 5, k — 10, and k = 15. We 
are stuck at k — 15, so we only need to show the cases k = 5 and k — 10. 

Let m be smallest such that [1, 2, 5, 5] does not represent m. Then 5 | m + 1 from 
above and |p^f is n °t a square modulo 5, taking £4 = 0. If r 7^ 0, then take £4 = 2, 
and note that 5 2r+1 (5n+ JO) — 2 • 3 • 5 = 20 (mod 2)5, so that m — 15 is represented by 
[1, 2, 5], and hence m is represented by [1, 2, 5, 5] as long as m > 15, which is as desired. 
Thus, we only need to consider r = 0. The non-squares modulo 5 are 2 and 3. As above, 
taking x 4 = 2 when 2m 5 +2 = 2 (mod 5) and 24 = 1 when 2 ™ +2 = 3 (mod 5) gives the 
desired result. 

For k — 10, we again conclude that r 7^ by taking x± — 1. Furthermore, if 2r ^ +2 = 3 
(mod 5), then x^ = 1 again gives us the desired conclusion. Hence, only the case 2m 5 +2 = 
3 (mod 5) remains. For [1, 2, 5, 10, k], with 5 | k, we may again conclude that by taking 
x 5 = and x 5 = 1 that 25 | k since 2(m + 2) = 5(5n + 2) and 2(m + 2 - k) = 5(5n' + 2), 
so 2k = 25(n — n'). Since 29 is the truant of [1, 2, 5, 10], the result follows when k 7^ 25. 
But we are stuck at k = 25, so we have the desired result. 
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We will leave out the analogous proofs for [1,2, 10], [1,3,4], and [1,4,6], but list the 
truants from their subtrees for completeness. The truants coming from [1,2,10] are 
29 and 49 (from [1,2,10,20]). The only truant from the subtree of [1,3,4] is 11. The 
truants from the [1,4,6] subtree are 17, 29 (from [1,4,6, 12], and 35 (from [1,4,6,6]). 

We will now show the subtree for [1,4,4]. While (1,4,4) is not genus 1, Benham, 
Earnest, Hsia, and Hung have shown that it is spinor genus 1 (0). Moreover they have 
shown that (1,4, 16) is spinor genus 1 by showing that the other member of its genus, 
namely Ax 2 + Ay 2 + 5z 2 + 4xz is spinor genus 1. Therefore, the difference r( lj4)4 )(8n + 
9) — r (i,4,i6)(8 ri + 9) can be decomposed into coefficients of the Siegel averaging of the 
genus, and a cusp form in U which has nonzero coefficients only at finitely many square 



classes. Using Schulze-Pillot's classification (1191 ) or the generalization of Earnest, Hsia, 
and Hung (0) to determine all t such that the square classe tl? has nonzero coefficients 
for the resulting cusp forms in U with these two quadratic forms, we conclude that only 
t — 1 occurs. Therefore, it follows that if m is represented by [1, 4, 4], then 8m + 9 must 
be a square. The first truant is m = 35, so we consider [1, 4, 4, k] for k < 35. Let m not 
represnted by [1, 4, 4, k] be given. Then m — k is also not represented by [1, 4, 4], so that 
8m + 9 is a square, say s 2 and 8m + 9 — 8k is a square, say t 2 . Therefore, s 2 — t 2 = 8k. 
But the difference between s 2 and t 2 must be at least the difference between s 2 and 
(s — l) 2 , which is 2s — 1. Therefore, 2s — 1 < 8A;. This restricts the possible choices for s 
to a (small) finite set, and hence the possible choices for m. Checking each such choice 
of s for each k allows us to determine the integers not represented by [1,4, 4, k]. We 
conclude that we are done for every integer other than k = 15 and k — 33. For k = 15, 
we represent exactly every integer other than 2 and 35, so [1, 4, 4, 15, will represent 
every integer except when k 1 = 33. For k — 33 or k' — 33, we are stuck at 35, and hence 
we are done with the [1,4,4] subtree. 

We will begin to use our GRH assumptions now. We will indicate clearly where we 
have made the assumptions, and furthermore we will indicate the cases where these 
assumptions were seemingly unavoidable. 

The eight cases [1,2,6], [1,2,8], [1,2,9], [1,2,11], [1,4,5], [1,4,8], [1,4,9], and [1,5,6] 
will all follow analogous arguments. The cases [1, 2, 6], [1, 2, 9], and [1, 4, 5] are the three 
cases where the GRH assumptions were seemingly unavoidable. In each case, we will 
be able to decompose the theta series for the corresponding quadratic forms with odd 
conditions into the Eisenstein series plus a Hecke eigenform. The Hecke eigenform in 
each case is in the complement of the space spanned by lifts of one dimensional theta 
series. Due to the fact that these are all genus 2 quadratic forms, we are able to first 
obtain the following proposition. 

Proposition 3.2. Each quadratic form [&i,&2>&3] in Tabled represents n if and only if 
the corresponding quadratic form 

ax 2 + by 2 + cz 2 + dxy + exz + fyz 

represents m. Here we do not have the "odd" condition for the quadratic forms. 
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Table 1. Equivalent Quadratic Forms 



Triangular Form 


represents n 


Quadratic Form 


represents m 


[1,2,6] 


/j 


(2,4,7,0,0,4) 


8n + 9 


[1,2,8] 


n 


(2,4,9,0,0,4) 


8n + 11 


[1,2,9] 


/( 


(2,3,4,2,0,2) 


2n + 3 


[1,2,11] 


n 


(1,6,8,0,0,4) 


4n + 7 


[1,4,5] 


n 


(1,4,5,0,0,0) 


8n + 10 


[1,4,8] 


n 


(4,4,9,0,4,0) 


8n + 13 


[1,4,9] 


n 


(1,4,9,0,0,0) 


8n + 14 


[1,5,6] 


n 


(3,3,4,0,2,2) 


2n + 3 



Proof. Each of the above assertions follows the same simple argument, which we will 
demonstrate explicitly for the form [1,2, 11]. 

First, note that if not all of x, y, and z are odd, then x 2 + 2y 2 + llz 2 = 8n + 14 has 
a solution modulo 8 only if x and z are both even. Therefore, if tb(n) is the number of 
solutions of the triangular form represented by b, and rqim) is the number of solutions 
to Q(x) = m, then 

£[i,2,n] 0) = ^(1,2,11,0,0,0) (8ra + 14) - r (4i 2,44,o,o,o)(8ri + 14) 

= ^(1,2,11,0,0,0) (8n + 14) - r (li2 ,22,o,o,o)(4ri + 7). 

Now we note that (1, 2, 11, 0, 0, 0) is a genus 2 quadratic form. The other representative 
of the genus is (2, 3, 4, 0, 02). If 2x 2 + 3y 2 + 4z 2 + 2yz = 8n + 14 has a solution, then it 
is clear that y must be even. Therefore, 

^(2,3,4,0,0,2) (8n + 14) = r (2,12,4,0,0,4) (8n + 14) = r ( i, 2 ,6,o,o,2)(4n + 7). 

But (1,2,6,0,0,2) is a genus 1 quadratic form, so it follows that r-(i j2i6j o,o,2)(4n + 7) is 
merely the value given by Siegel's local density formula as in (|23l). Using this observation 
and the fact that the local densities are equal for (1,2,6,0,0,2) and (2,3,4,0,0,2), it 
follows that r (lj2 ,6,o,o,2)(4n + 7) = r (2,3,4,0,0,2) (8rt + 14) = r (lj2 ,ii,o,o,o)(8n + 14). 

The form (1, 2, 22, 0, 0, 0) is again genus 2, and the other representative of the genus 
is (1, 6, 8, 0, 0, 4). Therefore, the theta series for Q = (1, 2, 22, 0, 0, 0), 

Oq(z) ■= Y.1 Qix) 

X 

satisfies 8q = Eq + g, where Eq is the Eisenstein series obtained by taking the local 
densities and g is a cusp form which is a Hecke eigenform. Siegel formula shows for 
Q' := (1,6,8,0,0,4) that 




1 



Therefore, 6qi = Eq — g. By obversing that the local densities for (1,2,11,0,0,0) 
and (4,2,44,0,0,0) are the same other than at p — 2, one can easily see by explicitly 
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computing the local density at p = 2 that aE (121100()) (8n + 14) = 2aE 1)2j83)0A0 (47i + 7). 
Therefore, we have shown that 

r(i,2,ii,o,o,o) (8n+ 14) -r (li2 ,22,o,o,o) (4n + 7) = 2a BlA22A00 (4n + 7) - (a ElA220M (4n + 7) +g) 

= a Ei,2, 22 ,o,o,o( 4ri + 7 ) - 9 = r (l,6,8,0,0,4)(4n + 7). 

This is precisely what we wanted to show. The other cases follow analogously. □ 

We now proceed by determining which integers in these arithmetic progressions are 
represented by these quadratic forms. Since each of these forms is genus 2, as well as 
spinor genus 2, we know that 9 = E + g, where g is a Hecke eigenform in the complement 
of the space spanned by lifts of one dimensional theta-series. Thus, we will employ the 
argument given in (fisl ) in the case where there is precisely one Hecke eigenform. 

We will begin by constructing an algorithm that given a non-zero Hecke eigenform g, 
a constant ce depending only on the local densities, such that the Eisenstein series has 
coefficients o,e{N) = CEh(mN) for some fixed integer m and N square free, where h(D) 
is the class number, the modulus q such that the corresponding twist from Waldspurger's 
theorem (I25I ) of the Shimura lift G of g is of modulus qN 2 , the integer m given above in 
the class number, an integer D such that a g (D ) is nonzero, a g (D ), and the character 
X that we are twisting by, and returning a bound D beyond which the Eisenstein series 
is nonzero in the congruence class corresponding to D . 

We do so by fixing X := 455, a : = 1.1573, and 02 '■= 1.3465 and calculating the 



bounds given in (Il5l ) for L(xmN, & + it) and L(G, %i a + it), where G is the Shimura lift 



of g. Bounding a, (3, 7, and S (these are all independent of the particular choice of g) 



as in Lemma 7.1 of (1151), we see that 



a(X) < 0.089028567932572, (3{X) > 0.0886630818642167, 
7(X) < 0.249235264918139, 6{X) < 0.0963544917482776. 

Given these bounds, the rest of the constants are easy to calculate (We use here the 
computer algebra system MAGMA (0).), and further details may be found in (fl3l ). 

Notice that the bound obtained by this algorithm is only valid in the congruence class 
congruent to Dq modulo 8 times the square of the determinant of the corresponding 
quadratic form, and that for each congruence class the choice of m and x ma Y vary. 
We therefore will run the algorithm for a choice of Dq in each congruence class which 
satisfies the congruence modulo a 2 power given above, and we will merely state the 
largest such bound obtained in Table [21 

We then check up to the bound D' for odd squares with a computer (using the fact 
that it is diagonal to our advantage by splitting off one dimension, as in (fl3l)) and list 
the integers not represented by each triangular form. We will list the triangular form, 
then the congruence classes not represented locally by the form, and finally the finite 
list of "sporadic" integers not represented globally but represented locally by the form. 
Note that we have a leaf (and hence GRH seems unavoidable with current techniques) 
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Table 2. Sufficient Bounds under GRH 



Triangular Form Bound Do for Quadratic Form Bound D' for odd squares 



[1,2,6] 


1.23 x 10 9 


1.23 x 10 9 


[1,2,8] 


6.0 x 10 s 


6.0 x 10 s 


[1,2,9] 


1.68 x 10 6 


6.72 x 10 6 


[1,2,11] 


8.0 x 10 4 


1.6 x 10 5 


[1,4,5] 


2.6 x 10 5 


2.6 x 10 5 


[1,4,8] 


5.7 x 10 s 


5.7 x 10 8 


[1,4,9] 


1.1 x 10 12 


1.1 x 10 12 


[1,5,6] 


4.00 X 1U 


1.82 X 1U 


Table 3. 


Exceptional Integers 


Not Represented 


Triangular Form 


Congruence Classes 


Sporadic Integers 


[1,2,6] 








[1,2,8] 





{4, 19,112} 


[1,2,9] 








[1,2,11] 





{4,25} 


[1,4, 5] 





{2,26,38} 


[1,4,8] 





{2, 16, 17} 


[1,4,9] 


m = 2 mod 9 and m = 8 


mod 9 


[1,5,6] 





{2,4,13,35} 



if and only if there are no congruence classes and no sporadic odd integers, namely the 
cases [1,2,6], [1,2,9], and [1,4,5]. 

Given table [3], it is easy to see that the only truants which arise from all of the subtrees 
other than [1, 4, 9] are 13, 17, 19, and 25. The first odd integer not represented by [1, 4, 9] 
is 11, so we only need to consider [1, 4, 9, k] for 9 < k < 11. The choice k = 9 is stuck at 
11, and k = 10 or k = 11 clearly represent every integer other than 2 and 8, by taking 
X4 = or z 4 = 1. 

The case [1, 2, 7, k] with 8 < k < 11 also follows from the above cases with no , and 
k = 7 is stuck at 11 so we get no new truants. The case [1,4,7] analogously gives no 
new truants, and [1,3,3] and [1,5,5] will follow similarly after we show the argument for 
[1, 3, 5] and [1, 5, 7]. Thus, it only remains to show the subtrees for [1, 3, 5] and [1, 5, 7]. 

For [1, 3, 5], we will follow a similar argument to above, but we must be slightly more 
careful. We note that 

t[i,3,5]( n ) = r (i,3,5,o,o,o)(8n + 9) - r ( i i5i i 2A o,o)(8n + 9). 

The theta-series for (1,3,5,0,0,0) decomposes as E + gi, while the theta-series for 
(1, 5, 12, 0, 0, 0) decomposes as \E + g 2 . In this case, the Shimuar lift of g\ — g 2 is 
\{— G15 + Ci 5 |\^(2) — 4Gi 5 |^(4)), where C7 15 is the (unique) newform of weight 2 and 
level 15 and V(d) is the <i-th V-operator (cf. (fl7l)). Therefore, since the Hecke operators 
commute with the Shimura lift, it follows that g\ — g 2 is a Hecke eigenform, so we may 
use the above algorithm with g — g\ — g%- 
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When 5 does not divide 8n+9, our algorithm gives the bound D' = 9.4x 10 s . Checking 
up to this bound we find no sporadic integers outside of n = 2 mod 5. Therefore, all 
integers other than a subset of n = 2 mod 5 are represented by [1, 3, 5]. However, since 
the truant of [1, 3, 5] is 7, we only need to consider [1, 3, 5, k] for 5 < k < 7, and 5 is 
stuck on 7. Hence, using i 4 = or i 4 = 1 we represent every integer other than 2. The 
reason for separating the case of n = 2 mod 5 is because the bound obtained by the 
algorithm was not computationally feasible and was not needed to obtain the desired 
result. 

In the case of [1, 5, 7], 

*[i,5,7](n) = ^(1,5,7,0,0,0) (8n + 13) - ^1,5,28,0,0,0) (8n + 13), 

and both (1, 5, 7, 0, 0, 0) and (1, 5, 28, 0, 0, 0) are genus 3 and spinor genus 3. Therefore, 
(1, 5, 7, 0, 0, 0) decomposes into 

E + gi + 02 

and (1, 5, 28, 0, 0, 0) decomposes into 

^E + g 3 + g 4 , 

where g\ , g2, g 3 , and g 4 are Hecke eigenforms in the complement of the space spanned by 
lifts of one dimensional theta-series. Moreover, computation shows that the Shimura lift 
of g 1 - g 3 is 1^E(G 35 + ciG 35 |^(2) + c 2 G 3b \V{A)), where G 35 is the newform of weight 



2V17 

2 and level 35 whose second coefficient is ~ 1 + v/ ^ . The constants C\ and c 2 are irrelevant, 
since we will twist the Hecke eigenform by killing all of the coefficients divisible by 2, 
and the coefficients of G\V(d) only has coefficients divisible by d. Moreover, the Shimura 
lift of g 2 - g A is (1 - l±ff){aG 3b + a( Cl )aG 3 ,\V(2) + a(c 2 )aG 35 \V(4)). Here a is the 

Galois map sending s/V7 to —\/Yf. 

We could write a separate algorithm from the one above for two (or any arbitrary 
number of) eigenforms, but for simplicity we will simply rewrite the above sum as 

(aE+( gi -g 3 ))+ ^-aj E+(g 2 

for some appropriate choice of a, and then bound each half separately as above, taking 
the maximum of the two bounds, since beyond the corresponding maximum bound Do, 
a aE+{ gi - gz ) > and i_ a \ E+ ^ g2 _ g ^ > 0. The optimal choice of a would give equal D 
bounds for each. 

We choose a — || here, so that | — a = ^. Doing so, whenever 5 does not divide 
8n + 13, we get a bound for the part corresponding to g\ — g 3 of D' = 4.53 x 10 8 and 
a bound of D = 1.65 x 10 8 for the part corresponding to g 2 — g 4 . The only sporadic 
integer not congruent to 4 mod 5 is 2. Therefore, we are again done, since [1,5,7, A;] 
represents every odd when 8<A;<9, asn = 4 mod 5 must be represented with x 4 = 1, 
and k = 7 is stuck at 9. 
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□ 

Remark 3.3. One may naturally ask here whether other currently known techniques 
are sufficient to determine the set So for S the set of all odd integers. This question, 
it turns out, is more easily asked than answered. We will give a partial answer to this 
question, however. For a node of depth 3, determining whether the node is a leaf when 
the node is not of genus 1 or spinor genus 1 does not seem possible given the current 
techniques. Indeed, if such an algorithm were known, it would most likely also give a 
solution for integers represented by Ramanujan's ternary quadratic form, and hence give 
an alternate proof of Ono and Soundararajan's result without GRH. 

There appear to be 3 leafs of depth 3 which are not (spinor) genus 1 when S is the 
set of all odd integers, as no counterexamples have been found up to the bound which 
is predicted by GRH. Therefore, if current techniques were sufficient to determine So, 
then these "leaves" must indeed not truly be leaves. But then there must exist (large) 
counterexamples, and our set So must be incorrect. For this reason, it is not this au- 
thor's belief that current techniques will solve this problem without GRH unless a new 
breakthrough occurs in the area. Moreover, since all of the other counterexamples for the 
other nodes (of greater depth) provably occur at very small integers, it seems reasonable 
to conjecture that So is as given. 

We conclude with the proof of Theorem 11.71 



Proof. The set S is precisely the integers represented by [2,3,4]. The integers less than 
or equal to 89 not represented by [2, 3, 4] are 1, 8, 31. Since the first truant of S is 2, we 
begin with [1] or [2]. 

The truant of [1] is 2, so we obtain [1,1] and [1,2]. For [1,2] the next truant is 4, 
and the remaining subtree has no further truants by Liouville's theorem. For [1, 1] the 
truant is 5. Each case other than [1, 1,3] is then covered by Liouville's theorem, while 
the case [1, 1, 3] has the same subtree as in Theorem 11.51 since in each case the truants 
were the smallest integer other than 8 not represented by the form and 31 is not one of 
the truants. From this we obtain the truants 17 and 89. This concludes the [1] subtree. 

For the [2] subtree the first truant is 3 so we escalate to [2, 2] or [2, 3]. The case [2, 2] 
still has truant 3, so we escalate to [2,2,2,3] and [2,2,3]. Clearly [2,2,2,3] represents 
every natural number greater than 1 because [2, 2, 2] represents every even natural num- 
ber. The form (2, 2, 3) is genus 1 so [2, 2, 3] represents every natural number n such that 
8n + 7 is not of the form 3 2r+1 (3m + 2). The first truant is 10 and if 3 { k, we are clearly 
done for [2, 2, 3, k], so we only have k — 3, 6, 9 and we are stuck at k — 9. For k = 3 we 
have the truant 19 and for k = 6 we have the truant 19. Analogous to the case [1, 1, 3] 
we are able to show that no other truants occur in this subtree. The truant of [2, 3] 
is 4, so we have [2,3,3] or [2,3,4]. We are stuck at k = 3 and [2,3,4] represents 5" by 
definition. □ 
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